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ABSTRACT 

We present a densely sampled early light curve of the optical/near-infrared 
(NIR) afterglow of the X-Ray Flash (XRF) 071031 at z=2.692. Simultaneous and 
continuous observations in seven photometric bands from g' to Ks with GROND 
at the 2.2 m MPI/ESO telescope on LaSilla were performed between 4 minutes 
and 7 hours after the burst. The light curve consists of 547 individual points 
which allows us to study the early evolution of the optical transient associated 
with XRF 071031 in great detail. The optical/NIR light curve is dominated by 
an early increase in brightness which can be attributed to the apparent onset of 
the forward shock emission. There are several bumps which are superimposed 
onto the overall rise and decay. Significant flaring is also visible in the Swift 
X-Ray Telescope (XRT) light curve from early to late times. The availability of 
high quahty, broadband data enables detailed studies of the connection between 
the X-ray and optical/NIR afterglow and its colour evolution during the first 
night post burst. We find evidence of spectral hardening in the optical bands 
contemporaneous with the emergence of the bumps from an underlying afterglow 
component. The bumps in the optical/NIR light curve can be associated with 
flares in the X-ray regime suggesting late central engine activity as the common 
origin. 
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Introduction 



Major progress in tlie understanding of the X-ray and optical afterglow light curves 
of Gamma-Ray Bursts (GRBs) and the softer X-Ra y Flashes (XRFs) has been made since 
the launch of the Swif t satellite (jGehrels et al.l 12004 ) and the rapid follow-up data provided 
by the X-Ray (XR T; burrows et all l2005bh and Ultra- Violet Optical Telescope (UVOT; 
Roming et al.l 120051 ) . However, access to the longer wavelength afterglow is still somewhat 
limited to the brighter h alf of all detected bu rsts. In contrast to the evidence of a generic X- 
ray afterglow light curve (INousek et al.ll2006l ). the few bursts with very e arly detected optica l 
counterparts s how considerable variety. For instance GRBs 990123 (lAkerlof et al.l Il999l ) 



and 041219A (iBlake et all l2005l : IVestrand et al.l l2005l : iMcBreen et all \200^ . have shown 



optical emission contemporaneous with the prompt phase of the burst. A significant delay 
of the apparent onset of the afterglow forward shock (FS) allowed ground based optical/NIR 



telesc opes to detect a risi ng component of the a fterglows for e.g. GRBs 03041 8 (IRykoff et al. 



2004\. 060418 060 607A flMolinari et al.ll2007l ). 070802 flKriihler et al.l l2008f ) and 071010A 
(ICovino et al.ll2008l ). A number of optical afterglows showed bumps superimposed onto the 
overall power law decay in late epochs, which are generally interprete d as the signature o f 
either inhomogeneities in the circumburst me dium (e.g. GRB 050502A.lGuidorzi et al.ll2005r) 
or late energy iniect ions (e.g. GRB 021004. Ide Uearte Postigo et al.lEoOsI or GRB 070311. 
Guidorzi et d1l2007h . 



The very early optical afterglow is of significant interest from a theoretical point of view 
and in particular how it relates to t he flares and plateaus seen in many X-ray afterglow 
light curves (e.g. O'Brien et al.l 2006 ). In the early phase the colour evolution is crucial to 
differentiate between different emission components. However, most of the rapid ground 
based follow up is obtained with robotic telescopes of small aperture size in white light 
or filter cycles. In both cases information about the spectral properties is absent or can 
only be obtained at relatively long times with respect to the dynamical time scale in the 
early evolution of GRB emission. The ambiguity between effects of a changing spectrum or a 
highly variable early light curve can only be addressed by systematic observations in different 
broad-band filters as synchronous and rapid as possible. Comprehensive da ta sets of early 
optical afterglows were published e.g. for GRB 021 004 Jlazzati et al.ll2002h . GRB 030329 



( Lipkin et al 



Racusin et al. 



2Q04\ . GRB 061126 flPerlev et al.ll2008f ) and the very bright GRB 080319B (e.g. 



20081 ) where the light curve is well sampled in time and frequency domains. 



suggesting that a standard jet break model alone can not account for the increasing variety 
of features in a GRB or XRF afterglow. 



Here we report on the optical follow up of GRB 071031 at redshift 2.692 (ILedoux et al. 



20071 ) using data obtained in seven broad-band filters from g' to Ks with the multi-channel 
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imager GROND ( Greiner et all2007l . l2008bl ). Ground-based optical/NIR observations started 
at ~4 minutes after trigger, yielding one of the best sampled early optical light curves. In 
combination with the detailed X-ray observations provided by the XRT, this constitutes a 
multi color light curve with spectral coverage from the near-infrared (NIR) to the 10 keV 
XRT band. 



2. Observations 



2.1. Swift 



The Burst Alert Telescope (BAT, iBarthelmy et al.l |2005| ) on-board the Swift satellite 
triggered on the long-soft G RB 071031 at To=01:06:36 UTC and immediately slewed to 
the burst (IStroh et al.ll2007bl ). The BAT light curve shows a two-peaked structure starting 
at Tq-IO s and ending at To+180s with a Tgo of 180±10 s. The fluence in the 15 keV to 
150 keV band is 9. Oil. 3 x 10 ~'' erg cm~^, with a flue nce ratio of 1.34 between the BAT 25- 
50 keV and 50-100 keV bands (jStamatikos et al.ll2007l ). This is remarkably soft compared to 
co nventional GRBs and q ualifies GRB 071031 as a XRF according to the working definition 
of ISakamoto et al.l (120081 ) . The BAT spectrum of the first peak is well described with a 
single power law with a photon index of 2.26±0.30 (x^=28.83 for 36 degrees of freedom). 
This is well outside the normal range for the low ener gy index a but si r nilar to the spectra l 
index f5 above the break energy for a Band function (IBand et al.lll993l : iPreece et al.ll2000l ). 
The peak energy of the prompt emission spec trum must then be close to or below the BAT 
lower energy range of around 30 keV (see also lMcGlynn et al.ll2005l : IStamatikos et al.ll2007l ). 
Therefore GRB 071031 is designated as XRF 071031 hereafter. 

The XRT began follow-up observations of the burst field 103 s after the trigger and 
detected an uncatalogued fading X-ray source at a position of RA(J2000)=00 ^ 25™ 37^.4, 
Decl(J 2000)=-58° 03' 33" with a refined 90% confidence error circle of 2'.'0 radius (jStroh et al. 
2007al ) . The early XRT light curve is dominated by significant amount of fiaring with bright 
fiares at around 120 s, 150 s, 200 s, 250 s and 450 s. Also the late X-ray data exhibit 
rebrightenings at 5.5 ks, 20 ks and 55 ks superimposed onto the overall power law decay. 
The complete XRT light curve is shown in Fig. [TJ 

The X-ray spectra were obtained with the xrtpipeline tool using the latest calibration 
frames from the Swift CALDB and standard parameters. The spectra were fitted with the 
XSPEC package (lArnaudl 119961) and a foreground hydrogen column density at the Galactic 
value of A^H=l-2xl020 cm"^ JXalberla et al.l 120051 ). 



The third instrument on-board Swift, UVOT started observations at To+114 s and 
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Fig. 1. — Swift XRT light c urve of the X-ray afterglow of XRF 071031 obtained from the 
XRT light curve repository (lEvans et al.l 120071 ) 
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found a transient source inside the XRT error circle in the white, v- and 6-band filters. The 
UVOT data show an incr ease in the brightness o f the afterglow of around 0.5 mag in the 
first few hundred seconds (IBreeveld &: Strohl 120071 ). 



2.2. GROND 



GROND responded to the Swift GRB alert and initiated automated observations that 
started at 01:10:21 UTC, 3 minutes 45 s after the burst and continued until local Sunrise at 
08:55:51 UTC. A predefined sequence of observations with successively increasing exposure 
times was executed and images were acquired in all seven photometric bands simultaneously. 
In total 84 individual frames in each g'r'i'z' and 1510 images of 10 s exposures in JHKs were 
obtained during the first night at airmasses between 1.1 and 2.4. The integration time 
of the CCD optical images scaled from 45 s to 360 s according to t he brightness of the 
optical afterglow. A variable point source was detected in all bands (IKriihler et al.l 120071 ) 
by the automated GROND pipeline (iKiipcii-Yoldas et al.ll2008l ) and its absolute position is 
measured to RA('J2 000)=00^ 25"] 37^ 24, Dec(J2000)=-58° 03' 33'.'6 compared to USNO-B 
reference field stars (jMonet et al. 2003) with an astrometric uncertainty of 0'.'3. Photometry 



and spectroscopy of the afterglow was also obtained by telescopes at CTIO (IHaislip et al. 



20071: ICobbll2007h and the VLT, the latter yiel ding an UVES and FORS spectroscopic redshift 
of 2.692 Jhedoux et al.l[2007l : IFox et al.ll2008h . 



Photometric calibration was performed relative to the magnitudes of ten secondary 
standards in the field of XRF 071031, shown in Fig. [2] and Tab. [TJ During photometric 
conditions, three spectropho tometric standard stars, SA114-750, SA114-656 and SA95-42, 
all primary Sloan standards (jSmith et al.ll2002l ). were observed with GROND. Observations 
of the GRB field followed within few minutes. The magnitudes of the Sloan standards were 
transformed to the G ROND filter system using their spectra and the GROND filter curves 
( IGreiner et al.l l2008bl ) . The obtained zeropoints were corrected for atmospheric extinction 
differences and used to calibrate the stars in the GRB field. An independent absolute 
calibration was obtained with respect to magnitudes of th e SDSS and 2MASS stars within 
the standard fields obtaine d from the SDSS dat a release 6 (lAdelman-McCarthy et al.ll2008l ) 
and the 2MASS catalogue (ISkrutskie et al.ll2006l ) with results consistent to the standard star 
calibration at the 0.03 mag level. 

Optical and i iear-infrared image reduction and photometry was performed using stan- 
dard IRAF tasks (lTodylll993l ). For each frame a model of the point spread function (PSF) 
was constructed using brighter field stars, and fitted to the afterglow. The relatively large 
seeing between 2" and 3" together with the pixel scale of 0'.'16 for g'r'i' z' and 0'.'60 for 
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Fig. 2. — GROND r' band image showing the afterglow of XRF 071031 and the XRT error 
circle. The secondary standards used for calibration are labeled 1 to 10 and listed in Tab HJ 
The lower left corner shows a zoom in to the afterglow position. 
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JHKs resulted in an excellent spatial sampling of the PSF with statistical fit errors of order 
0.2 % for g'r'i' z' and 0.5 % for JHKs ■ For consistency, we also performed standard aperture 
photometry with compatible results with respect to the reported PSF photometry. All data 
were corre cted for a Galactic f oreground reddening of Eb-v = 0.012 mag in the direction of 
the burst JSchleeel et al.lll998h . 

The stacking of individual images was done twice for different purposes. Firstly, all 
available data were used and individual frames were stacked until a statistical error in the 
PSF fit of around 0.1 mag was obtained. This resulted in 75 frames in each g'r'i' , 51 in 
z' , 118 in J and 76 in each H and Kg, yielding the multi-wavelength light curve shown 
in Fig. [31 Secondly, only NIR data simultaneous to the optical integrations were selected, 
excluding the frames which were taken during the g'r'i' z' CCD read out. The resulting NIR 
image stacks were used to derive the optical to NIR spectral energy distribution (SED). 



3. Analysis 

3.1. The optical/NIR light curve 

The multicolour light curve of XRF 071031 is complex and not described by smoothly 
connected power laws alone (Fig. [3]). Evident in all colors is the initial increase in brightness 
up to ~1 ks, which smoothly turns over to a generic power law decline with superimposed 
variations. In order to study small scale irregularities at the highest possible signal-to-noise, 
a white-light curve was derived by combining all g'r'i' z' CCD data (Fig. HI upper panel), 
which are exactly synchronous by hardware setup. The NIR bands are excluded from this 
process due to the intrinsically larger photometric error which would serve only to increase 
the uncertainties in the summed data. 

In order to better visualize the afterglow light curve, it is also presented differentiated 
in its native log-log scale, i.e. S{log{F^))/6(\og(t)), which directly represents the local power 
law decay index a(t) (Fig. HI lower panel). On a log-log scale, the first derivative of a power 
law is a constant, so for a smoothly connected power law rise and decay, one would expect a 
positive constant at early, turning into a negative at later times. Changes in the power law 
index and deviations from the decay are clearly visible in this representation. 

Without imposing an a priori model of the afterglow. Fig. H] convincingly demonstrates 
that the overall trend of the light curve i s well d escribed by two smoothly connected power 



laws as introduced in iBeuermann et al.l (119991 ). In both panels of Fig. HI however, the 



deviations from a Beuermann-like power law rise and decay are clearly apparent. There are 
two features, which require either a superimposed component, or a different parametrization 



- 8 - 




Fig. 3. — GROND optical and NIR light curve of the afterglow of XRF 071031 taken between 
~4 min and 7 hours after the trigger. 
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of the intrinsic afterglow. The first and most obvious is the additional emission component 
in regions A,B,C,D shown in the upper panel of Fig. |H This extra emission requires features 
intrinsic to the source or its environment to produce the observed flux excess with respect to 
the power law. The second is the steepening of the power law decline after a rebrightening 
from region I over II to III in the lower panel of Fig. HI However, there is so much variability 
within the light curve that the underlying afterglow cannot be established with high certainty. 

The early rising component in the optical bands might be related to the deceleration of 
the FS by the circumburst medium, which happens when the swept up medium efficiently 
decelerate the ejecta. From the time of the light curve peak, the initia l bulk Lorent z facto r 
of the outflow Fq can be constrained. Using the formalism outlined by ISari &: Piraru (119991 ) 

Panaitescu fc Kumail J200oh and iMolinari et aP J2007h . Fq is estimated to ^ 90 



in a wind shaped circumburst medium and to ~ 200 



E53 

rio.2n 



1/8 



'70.2^* 



1/4 



for an ISM type environment. 



with a weak dependence on the uncertain parameters A* being the normalized wind density, 
£^53 the isotropic-equivalent energy released in 7-rays in 10^^ erg, 779.2 the radiative transfer 
efficiency normalized to 0.2 and n the ISM density in cm~^. The slow rise with a power law 
index of ~0.7 suggests a wind like environme nt, which would on ly be consistent with the 
closure relations for a very hard electron index (jPai &: Chengll200ll ) of p~1.6 and < v using 
the spectral and temporal slopes i n the late afterglow light curve ax,o ~1 and /?x=0.8±0.1. 



The classical closure relations (e.g. IZhang &: Meszarosll2004l ). however, would favour an ISM 



environment in the slow cooling case with and a more canonical value of p~2.6. 



[e.g. 



Alternatively, the initia l rise could be the result of a structured outflow seen off-axis 

I n the case of previous fast and slowly rising afterglow light 



Panaitescu et al. 



curves, 



199 



Panaitescu fc Vestrandl (120081 ) flnd an anticorrelation of peak flux in the R band F^^r 
and peak time tp. K-correcting the afterglow to 2; = 2 to match the previous sample, we 
flnd that the optical/NIR light curve of XRF 071031 flts very well into this anticorrelation. 
In this interpretation the slow rise would hint on a smooth angular structure of the outflow 
( jPanaitescu fc Vestrandl l2008l ). 



Chromaticity around peak brightness was tested by comparing the optical/NIR SED 
before and after the total maximum. Apart from changes in the spectral index which can 
be attributed to the emergence of the bumps (see Sec. 13.41 and Fig. [7]), there is no evidence 
for a change in the spectrum before and after the light curve peak. The time of the light 
curve maximum is not correlated with energy and all bands peak at a similar time within the 
measurement uncertainties. Such evolution would be expected if the main peak was caused 
by cooling of the ejecta after the prompt emission res ulting in the shift of th e characteristic 
synchrotron frequency into the optical bands (e.g. IZiaeepour et al.l 120081 ). In addition, a 
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moving Urn through the opti cal bands is expected to cause a stron g change from a positive to 
negative spectral index (e.g. lSari et al.lll998l : lGranot fc Sarill2002l ). Neither effect is observed. 



After correcting for Galactic foreground extinction, the SED from g' to Ks was fit 
with a power law and dust reddening templates in the host at 2;=2.692. Using extinction 
templates from the Milky Way (MW), Large- (LMC) and Small Magellanic C loud (SMC) and 
supernovae induced dust (Pej 1992 : Maiolino et al.lEo04 : Stratta et al. 2007 ) no signatures of 
dust extinction in the GRB host are evident down to la confidence limits of A^'^^ < 0.06 mag 
(MW), A^"""^ < 0.07 mag (LMC), A^"""^ < 0.05 mag (SMC and SNe induced). Therefore the 
effect of dust reddening is considered as negligible in the following analysis. The deviation 
from a power law SED in the GROND g' band is consistent with Lyman-a absorption in the 
GRB host at ^=2.692. 



3.2. The X-ray afterglow light curve 



In addition to the variable and densely sampled light curve in the GROND filter bands, 
the X-ray afterglow is bright and well covered by XRT observations. Similar to what is seen 
in the optical bands, the X-ray data show strong variability, and the underlying afterglow 
is poorly constrained. After excluding the very early data, where there are no GROND 
observations (t-To < 300 s), we fitted the remaining data using a similar procedure as used 
for the optical bands with a combination of a smoothly connected power laws. Remarkably, 
the obtained late power law index ax=0.99±0.12 is compatible with the best fit from the 
GROND data ao=0.97±0.06, providing additional evidence that the applied model fitting 
traces the underlying power law decay of the afterglow reasonably well. The XRT light 
curve and the residuals to the power law fits are shown in Fig. O After flaring episodes, the 
XRT light curve drops back to the power law, consistent with the rebright enings observed 



in previous GRBs or XRFs (e.g. iBurrows et al.ll2005al : iRomano et al.ll2006bl ). A H XRT light 



curve data have been obtained from the Swift XRT light curve online repository fevans et al. 
2007h . 



3.3. The bumps 



Both X-ray and optical/NlR data show significant variations from the typical power 
law profiles at very early times. While this is observed in nearly 50% of all X-ray afterglows 
Nousek et al.ll2006l ). this is rarely seen so early in an optical light curve starting ~400 s 



le.g. 



after the trigger. 
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Fig. 4. — White light curve of the afterglow of XRF 071031 in the upper panel. A numerical 
differentiation of the data in the native log-log scale is shown in the lower panel. Regions of 
bumpiness are denoted with capital letters (upper panel) , and the steeping of the power law 
in Roman numbers (lower panel). The effect of a changing smoothness parameter s of the 
turnover is indicated by different lines in both panels. 
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Fig. 5. — XRT light curve of the afterglow of XRF 071031 in the time domain simultaneous 
to the GROND data (top panel). In the middle panel the residuals to the smoothly connected 
power law fit, as well as their modehng using Gaussians. The residuals of the first flare have 
been scaled by 1/500. The bottom panel shows the hardness ratio of the 1.5-10 keV versus 
0.3-1.5 keV band. 
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We fitted the combined white light data with a canonical afterglow rise and decay with 
superimposed Gaussian profiles to account for the variations in the light curve. While this 
may not be the true physical model, it provides a good fit to the data, and represents the 
morphology of the bumps with adequate accuracy. For the very early optical data, where 
the time sampling of the light curve is naturally sparse, we used a cubic spline interpolation 
with equally spaced nodes in the native log- log scale of the afterglow to constrain the fit. 
The data suggest the existence of three major and three small bumps: The three brighter 
ones, TTi to 773, which peak at 0.6 ks, 6 ks and 18 ks, and three fainter ones, ^1 to ^3, all shown 
in the lower panel of Fig. [61 The first faint bump is only indicated by one data point and 
the last peak is not sampled by the observations due to the break of dawn. The fit shown 
in Fig. [6] has a Xred ~ ^-^^ with 52 degrees of freedom (dof). The best fit parameters are 
presented in Tab. [21 

The typical timescale of variation for the optical bumps is (AT/T) = 0.35 ± 0.13, where 
AT is calculated as the full width half maximum of the Gaussian, and have a relative flux 
increase with respect to the underlying afterglow of (AF/F) = 0.29 ± 0.18. All bumps, 
which are reasonably sampled by our observations, have a change in the slope 6a between 
0.5 and 0.9, which can be directly deduced from the lower panel of Fig. [H 

Due to the faintness of the source, the late X-ray afterglow light curve is not equally 
well sampled as the optical and a detailed analysis is not possible in late epochs. There is, 
however, evidence for three X-ray flares (pi to ^3 between 350 and 25000 s with synchronous 
optical coverage. To derive the morphology of the flares in the XRT light curve, we use a 
similar approach as for the optical bands for a direct comparison: An underlying continuum 
fitted with power laws and superimposed Gaussians to account for the evident flux excess 
which yields a reduced xled of 1.15 for 154 degrees of freedom. The best fit is shown in Fig. [51 
and the corresponding parameters are reported in Tab. [21 

The X-ray flares are much stronger (AF/F) = 1.28 ± 0.28 with respect to the underlying 
afterglow t han the optical bump s. Comparing against a statistical sample of previous X- 



ray flares (IChincarini et al.l 120071 ) shows that the flares observed in XRF 071031 populate 
a similar phasespace region of AT/T versus AF/F and thus resemble the morphology of 
previous flares. 



3.4. Spectral evolution 

The spectral index j3, where F^, oc u"^ of the optical/NIR SED is observed to evolve 
with time (Fig. [Tj). Remarkably, the chromatic evolution is correlated with the residuals 
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Fig. 6. — White liglit curve of tiie afterglow of XRF 071031 (upper panel, shown are the data 
and the cubic spline interpolation with equally spaced nodes in log(t)). The data were fitted 
using the sum of a smoothly connected power law for the canonical afterglow (dashed line) 
and Gaussian profiles to account for the evident flux excess (solid line). In the lower panel 
the residuals to the smoothly connected power law, as well their modeling by six Gaussians. 
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of the data against the hght curve fits. Applying a standard statistical correlation analysis 
yields a correlation coefficient of ~-0.73, and thus a null hypothesis probability of ~ 10~®. 

The correlation of spectral hardening and bumpiness suggests, that this is an intrinsic 
feature of the emission component in the bumps, rather than the afterglow itself. In the 
ISM model for example, the cooling fr equency z^^ ; movi ng through the optical bands, would 
identify itself by a spectral softening (ISari et al.l Il998l ) , in contrast to the o bservations. A 
hardening of the spectrum would be expected in a wind like environment (IGranot fc Sari 
20021 ) with a change in the spectral index of 0.5. The observed change (see Fig. [Tj) is ~0.2 
and therefore not compatible with a cooling break passing the optical bands also in the wind 
model. 

In addition to the observed chromatic evolution in the optical bands, also the X-ray 
data show strong spectral changes throughout the observations. This is already indicated 
by the evolving hardness ratio of the two XRT bands (1.5-10 keV over 0.3-1.5 keV), shown 
in the lowest panel of Fig. [51 

This evolution becomes more evident, when including individual early flares from Fig. [T] 
into the analysis. For each flare in Tab. HI a spectrum was extracted and fltted by single 
and broken power law models. The rest frame column density A^h at 2=2.692 was obtained 
by fltting the late photon counting data and found to be consistent with zero within a la 
confldence level. Combined with the negligible Ay from the optical data and assuming a 
constant Ah in the burst environment, the intrinsic Ah is neglected in the spectral flts. There 
is strong evidence that the flares are better modeled by broken rather tha n single pow'er law s 
as shown in Ta b. HI which is simila r to flares see n in e.g. GRB 05 1117A (IGoad et al.l 120071 ). 
GRB 050713A dGuetta et al]l2007h and 061121 JPage et al.ll2"o07l l 



From the X-ray data alone, there is evidence for a break in the spectrum in th e 1 keV 
range for the early flares. This is consistent with the result of iButler fc Kocevskil (120071 ). 
who flnd that the peak energy of the flare spectrum Ep crosses the X-ray bands on a tj^ical 
timescale of 10^ to lO'^ s. Combining the excess flux in the optical bands for tti and the 
X-ray spectrum in 0i, the spectrum can be constrained over a broad energy range. For 
the optical bands the dominant emission process is FS emission even at early times, and 
the afterglow model fltting was used to disentangle the different components. In this way 
estimates of the ffux attributed to the ffare component can be obtained. As shown in Fig. [HI 
the broadband spect rum of the first ffa re is reasonably well (x^=129 with 114 d.o.f) described 
by a Band function ([Band et al.lll993l ) with an Ep of 1.79 ± 0.59 keV and a very reasonable 
set of parameters q ;=— 0.78i:0.03 and /9=— 1.921;q;}7 as compared to the BATSE sample 
jKaneko et al.ll2006h . 
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Fig. 7. — Color evolution of the optical/NIR afterglow of XRF 071031. A direct comparison 
of the spectral index (5 (upper points - left y axis) against the residuals to a fit of a canon- 
ical afterglow light curve model (lower points - right y axis). The plot clearly indicates a 
correlation of the spectral hardening and the rebrightening episodes in the light curve. 
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We caution that this fit imphcitly assumes that the excess emission seen in the optical 
bands is correlated with the X-ray flare and the applied model of the underlying FS emission 
traces the afterglow reasonably well. As the fundamental shape of the afterglow can be 
different then the empirical Beuermann profile, this might introduce significant systematic 
errors in the analysis. Additionally, an underestimated or even varying column density A^h, 
would change the soft X-ray absorption and thus the overall broadband and X-ray fits. 

For the later flares the different components attributed to FS and flare emission cannot 
be deconvolved with high certainty. The X-ray data are faint and affected by the underlying 
afterglow, and the excess flux in the optical bands is strongly dependent on the parameters of 
the light curve fitting. In particular, there is a strong ambiguity between the sharpness of the 
break and the light curve decay. Given the large uncertainties, the optical colour attributed 
to later bumps is compatible with the slope of the lower energy part of the Band function 
found for the first flare. Ep, however, can no longer be constrained, but it is interesting 
to speculate that a similar Band function with an Ep between the optical and X-ray bands 
would account for the observed excess fluxes. 



4. Discussion 



4.1. The likely cause of the flares 

Previously, bumps or flares in optical afterglow light curves have been reproduced us- 
ing either a superimposed reverse shock comp onent for early flares, inhomogeneities in the 
circumburst medium (e.g. IWang fc Loebl 120001) or the angular distribution of the energy in 
the jet (pat chy shell model, e.g. iKumar fc Piranll2000l ) or late energy injection by refreshed 



shocks (e.g. Rees fc Meszaros 19981 ) for later flares. However, a clear discrimination in the 



individual previous cases was not possible due to the lack of simultaneous coverage in broad 
wavelength ranges. 

Although we cannot completely rule out a reverse shock component for the first bump, it 
seems likely that it is produced by the same mechanism as the later ones. After subtracting 
the rising power law, the decline of the first bump can be fitted with a power law of index 
~-7. This would be surprisingly fast f or a reverse shock, wh ich is expected to dechne with 



t ^ in the basic fireball model (e.g . iNakar fc PiranI |2004| ) and not faster than ~ t for 



more complicated models (IKobayashil |2000| ) . The fact that the morphology is comparable 
to the later bumps additionally hints on a common origin. Therefore we try to account 
for all bumps with as few assumptions as possible, thus searching for a phenomenological 
explanation for all bumps observed. 
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Fig. 8. — Broadband spectrum of the first fiare tti with simultaneous optical and X-ray 
coverage. The data are fitted by a Band function with Ep=1.79±0.59 keV, Q!=— 0.78±0.03 
and p=-1.92l:lfT 
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Recently, iNakar fc GranotI (120071 ) found, that density jumps in the circumburst medium 
cannot account for the majority of fluctuations previously observed in GRB afterglow light 
curves. Only large contrasts in the circumburst medium density are able to produce bumps 
with a change in the temporal power law decay index 5a ~1 in the light curve with a long 
transition time, which scales linearly with contrast. Thus, signatures of inhomogeneities in 
the circumburst medium in the optical light curve are expected to be smoother to what is 
observed. 

A jet with inhomogeneities in the angular energy distribution produces episodic bumps 
when the coiie of a relativistically beamed patch enters the field of the view of the observer 
( iZhang et al.ll2006l ). Similar to refreshed shocks, these patches inject additional energy into 
the blastwave, the afterglow er aission is boosted to a high er level and resumes the same power 
law index as before the bump (jZhang fc Meszarosll2002l ). A characteristic for rebrightenings 
due to patchy shells or refreshed shocks consequentially is a step like afterglow light curve. 
Given the steepening of the power law post bump (Fig. Hj) this scenario seems inconsistent 
with the bumps observed in XRF 071031. We note however, that the refreshed shock scenario 
can produce rebrightenings on relatively short ti mescales under certa in conditions, where the 
light curve drops back to the initial decay (e.g. iGuidorzi et al.l 120071 ). 



Although the morphology of the light curve and optical bumps in the GROND data 
might be explained within the framework of variable external density or energy dissipation 
in the FS, both the spectral evolution as well as the correlation with the X-ray data argue for 
an independent origin and a second emission component . A hardening of the spectral index 
for flares in an optical light curve is not unprecedented (iGreiner et al.ll2008al ). and suggests 
a different emission than the generic afterglow FS. In addition, the first pronounced bump 
in the optical light curve is already observed during the rise of the afterglow at 600 s, when 
the apparent onset of the FS just started. 

If all previously observed X-ray fiares are due to the same physical process, it is 
very likely that they do not originate from external s hocks that giye rise to the after- 



glow emission, but from late time internal shocks (e .g. iBurrows et al.l l2005al : IZhang et al. 



2006 



Chincarini et al.l 120071 : iButler &: Kocevskil |2007| ). Thus, X-ray flares seem to be pro- 
duced by a similar mechanis m as the prompt 7-ray s, which are also caused by internal 



shocks in the standar d model (IRees fc Meszaroslll992l ). Detailed analysis for GRB 050820 A 
(IVestrand et al.ll2006l ) showed that the optical emission contemporaneous with the prompt 
phase can be explained as the superposition of forward shock and emission correlated with 
the 7-rays. Similarly for XRF 071031, after subtracting the dominating FS in the optical 
bands, the flare spectrum from NIR to X-rays is well described with a Band function. 



Remarkably, the optical bumps show features which have been previously observed in 
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X-ray flares: a hardening of their spec tra (e.g. lBurrows et al.ll2005al : lButler fc Kocevskill2007 
Falcone et al.l 120071 : iGoad et al.ll2007l ) and a correlati on of the duration with the time where 
the b ump occurs, i.e. a roughly constant AT/T (jChincarini et al.l 120071 : iKocevski et al. 



20071 ). They are however, less pronounced than typical X-ray flares. Temporal analysis of 
the early data shows, that the optical bump peaks significantly later than the X-ray flare. A 
hard to soft evolution of Ep whi ch is found in the majority of all bright flares where a detailed - 
spectral analysis is p ossible (e.g. Burrows et al.|[2005a : Romano et al. 2006a : Perri et al. 2007; 
Falcone et al.l 120071 ) provides a natural explanation for the time difference between the flare 
in the X-ray and optical wavelength range. Spectral lags and a broadening towards lower 
ener gies have been observed in a number of previous X-ray fla res and the prompt emission 



[e.g. 



Norris et al.lll996l : iRomano et al.ll2006bl : IPerri et al.l 120071 ) 



Based on the light curve fitting, the observed peak of the early optical flare is delayed 
by T ~ 130 s compared to the X-rays, which corresponds to 35 s in the bursts rest frame. As 
the temporal coverage of the optical light curve is sparse in the early time frame, this delay 
is strongly dependent on the assumed functional form of the afterglow and flare morphology, 
but in any case it is significantly longer than what is typically observed as spectral lags 
in the prompt phase. Typical values for prompt lags range from sli ghtly negative (i.e . 



soft preceding hard bands) to several seconds for long la g GRBs (e.g. iNorris et al.l 12000 
Gehrels et al.ll2006l : iMcBreen et al.ll2008l : iFoley et al.ll2008l ). The observed time difference in 
the case of XRF 071031, however, is based on entirely different energy ranges. In particular, 
the difference between hard and soft energy bands is around a factor of 10^ for X-rays versus 
optical, while it is ~10 for BAT channels. 

Combining bump morphology, colour evolution, broadband spectrum and the temporal 
connection to the X-ray data as shown in Fig. [9], the most likely origin of the bumps in the 
optical/NIR light curve of XRF 071031 is the same as in the X-rays, namely the soft tail of 
emission correlated with late internal shocks. 



4.2. Paucity of detection of correlated early optical bumps and X-ray flares 



While a nur nber of bursts have shown considerab l e variability in the optical bands durin g 
X-ray flares (e.g. IPandev et allbood : iBoer et al.l[2006l : IPage et al.l[2007l : iMalesani et"aDl2007h , 
the early optical light curve is generally not as strongly affected as the X-rays by flaring 
episodes. The UVOT and ground based robotic telescopes with their fast response provide 
a good sample of optical observations simultaneous with the XRT light curve. If the first 
optical bump and the first X-ray flare in XRF 071031 are correlated as suggested by the 
data, the signature of a typical X-ray flare in the optical light curve can be approximated 
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1000 10000 



Time since trigger [s] 

Fig. 9. — Combined plot of XRT light curve (crosses, left y axis) and GROND white light 
(dots, right y axis). The vertical dashed lines indicate the emergence of the bumps from the 
underlying optical afterglow emission. 
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with the ratio of absolute optical to X-ray peak flux obtained for XRF 071031. 



GRB 060607A for example, showed significant X-r ay flaring very early in the X-ray light 



curve with only marginally correlated optical emission (iZiaeepour et al.ll2008f). However, the 
afterglow was very bright, reaching ~14.3 mag in the r band (INysewander et al.l 120071 ). If 



the flare was described with a similar spectrum as in XRF 071031, the expected optical 
counterpart of the brightest X-ray flare is estimated to have a peak amplitude of ~200 /iJy 
and could only be resolved with continuous photometric monitoring with a precision of at 
least 0.03 mag. Thus, if the emission in the flares is not strong enough with respect to the 
FS, a bright afterglow can easily outshine flare signatures in the optical bands even for very 
bright X-ray flares. In fact, GRB 061121 (jPage et al.ll2007l ). had simultaneous coverage with 
BAT, XRT and UVOT of the prompt emission with a peak flux density in the X-rays of ~15 
mJy. iPage et al.l (120071 ) find, that most of the flux is emitted in the 7-ray and X-ray bands, 
while the UVOT data only showed a relatively small increase in brightness. 



5. Conclusions 

The detailed light curve and spectral evolution of the afterglow of XRF 071031 showed 
several remarkable features. A general trend of an achromatically rising and falling opti- 
cal/NIR afterglow, consistent with a number of previous bursts where an early increase in 
brightness was reported. The achromatic turnover from rise to decay suggests the apparent 
onset of the FS as its origin, either due to the increase in radiating electrons in the pre- 
deceleration phase or a structured outflow seen off axis. In the latter case the increase in 
brightness is caused by the deceleration of the shock, so that the relativistically beamed cone 
widens and gradually enters the sight line of the observer. 

Superimposed onto the afterglow continuum emission are bumps which have a harder 
SED and appear in similar epochs in the optical/NIR as well as in the X-ray bands. Based on 
the GROND data collected for XRF 071031 we conclude that the variation in the light curve 
are most probably the optical counterparts of X-ray flares, and therefore directly connect 
variability in an optical light curve with the emergence of X-ray flares. Although emission 
from external shocks or a combination of different effects cannot be completely ruled out, an 
internal origin seems to nicely account for the majority of observations: The light curve shape 
and in particular the morphology of the bumps, the spectral hardening in the optical SED, 
the observed decrease in E-p from the prompt emission to the flares and the broadband flare 
spectrum from NIR to X-rays. The spectral similarities to the prompt phase strengthen the 
picture of X-rays flares as later and softer examples of the prompt emission due to internal 
shocks. This connection provides additional evidence that inner engine activity may last or 
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be revived on a timescale of hours or days at least for some bursts. 
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Table 1. Secondary standards in the GRB field in the GROND filter bands from 1 to 10 



Ra/Dec 


g' 


r' 


i' 


z' 


J 


H 


Ks 


Deg [J2000] 


mag 


mag 


mag 


mag 


mag 


mag 


mag 



6.4098/-58.0831 15.44 ± 0.02 14.55 ± 0.02 14.28 ± 0.02 14.13 ± 0.02 13.09 ± 0.04 12.55 ± 0.04 12.50 ± 0.06 

6.4186/-58. 072.3 16.32 ± 0.02 15.73 ± 0.02 15.54 ± 0.02 15.43 ± 0.02 14.46 ± 0.04 14.07 ± 0.04 14.00 ± 0.06 

6.3933/-58.0662 17.84 ± 0.03 17.32 ± 0.03 17.18 ± 0.03 17.10 ± 0.03 16.18 ± 0.04 15.94 ± 0.06 15.65 ± 0.07 

6.3672/-58.0536 20.71 ± 0.05 19.33 ± 0.03 18.53 ± 0.03 18.13 ± 0.03 16.86 ± 0.05 16.41 ± 0.07 16.28 ± 0.08 

6.4579/-58.0492 20.03 ± 0.04 18.90 ± 0.03 18.49 ± 0.03 18.31 ± 0.03 17.07 ± 0.05 16.68 ± 0.07 16.45 ± 0.08 

6.3663/-58.0782 20.09 ± 0.04 18.62 ± 0.03 17.60 ± 0.03 17.13 ± 0.03 15.80 ± 0.04 15.30 ± 0.05 14.92 ± 0.06 

6.3751/-58.0631 21.26 ± 0.05 20.13 ± 0.04 19.72 ± 0.04 19.52 ± 0.04 18.29 ± 0.07 17.87 ± 0.08 17.79 ± 0.10 

6.4021/-58.0496 22.72 ± 0.10 20.85 ± 0.05 20.06 ± 0.05 19.65 ± 0.05 17.94 ± 0.06 17.35 ± 0.07 16.11 ± 0.08 

6.4117/-58.0575 22.49 ± 0.08 21.21 ± 0.06 20.76 ± 0.06 20.51 ± 0.07 18.56 ± 0.07 17.78 ± 0.08 16.86 ± 0.10 

6.4375/-58.0485 21.28 ± 0.05 20.77 ± 0.05 20.61 ili 0.05 20.59 ili 0.07 19.29 ± 0.09 18.91 ± 0.10 17.98 ili 0.15 



Table 2. Parametrization of the excess flux in the GROND bands using Gaussians 



Bump 




AT/2 [s] 


Amplitude [/iJy] 


AT/T^,d 


AF/F 


TTl 


587 ± 12 


84 dz 18 


149 ± 5 


0.29 ± 0.06 


0.24 ± 0.01 




1726 ± 6 


130 ± 4 


39 ± 5 


0.15 ± 0.01 


0.06 ± 0.01 


7r2 


6302 ± 156 


1648 ± 136 


60 ± 5 


0.52 ± 0.05 


0.26 ± 0.02 




9439 ± 377 


1461 ± 254 


29 ± 4 


0.31 ± 0.06 


0.18 ± 0.02 




16505 ± 135 


4013 ± 263 


43 ± 2 


0.49 ± 0.03 


0.46 ± 0.02 




28000 


5000 


32 ± 4 


0.36 


0.55 ± 0.03 



Table 3. Parametrization of the excess flux in the X-ray bands using Gaussians 



Bump 




AT/2 [s] 


Amplitude [Counts/s] 


AT/T^id 


AF/F 




459 ± 2 


42 ± 2 


14.5 ± 0.8 


0.18 ± 0.01 


1.36 ± 0.07 


</)2 


5528± 85 


604 ± 114 


0.047 ± 0.01 


0.22 ± 0.04 


1.20 ± 0.20 


4>3 


20507 ± 952 


4557 ± 1088 


0.012 ± 0.002 


0.44 ± 0.10 


1.28 ± 0.13 



Table 4. Spectral fits to the XRT flares using XSPEC. 



Flare Times [s] Model Photon Index 1 Break energy [keV] Photon Index 2 /doi 



1 
1 


137 - 
137 - 


150 
150 


Single 
Double 


1 r,o + 0.05 

^^■^■^-0.05 
r, 79 + 0.09 
'^-0.09 




1 85+°-l2 


323/131 
157/129 


2 
2 


189 - 
189 - 


217 
217 


Single 
Double 


-, -9+U.U9 




2 20+° " 


180/85 
89/83 


3 
3 


236 - 
236 - 


284 
284 


Single 
Double 


1 07+U.U5 

1 90+0.19 
^^■^■^-0.18 


J 07+0.16 


9 oq+0.17 


146/74 
67/72 


4 (*i) 
4 (*i) 


396 - 
396 - 


547 
547 


Single 
Double 


1 64+"'"* 

1 26+°-^^ 
^■■^"-0.32 


0.98t°;^g 


1 RS+008 
i-»o_0.08 


172/110 
120/108 



